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Abstract We investigate the pseudo-spin symmetry case of a spin- 1 particle governed by 
the generalized isotonic oscillator, by presenting quasi exact polynomial solutions to Dirac 
equation with pseudospin symmetry vector and scalar potentials. The resulting equation 
is found to be quasi-exactly solvable owing to the existence of a hidden si (2) algebraic 
structure. A systematic and closed form solution to the basic equation is obtained using the 
Bethe ansatz method. Analytic expression for the energy is obtain and the wavefunction 
■j , is derived in terms of the roots to a set of Bethe ansatz equations. 

-5 



PACS 03.65.-w, 03.65.Fd, 03.65.Pm, 03.65.Ge 

7-H \ 

Keywords: Generalized isotonic oscillator, Bethe ansatz method, Quasi-exactly solvable 
■^j- ■ systems, Dirac equation 

(N 

in 

^ ■ 1 Introduction 

O 

CN ' In a recent study [I], the nonrelativistic one-dimensional quantum system described by 

the nonpolynomial potential of the form 

w 2 x 2 x 2 — a? 

X: v{x) = — + 9a WT^f (1) 

where a is a positive parameter, was shown to be exactly solvable in the particular case 



of g a = 2 and wo =1/2, with a general solution 

Pu(x) 

(2x 2 + l) 



^ ( x ) - e~Z 2 /2 



(2) 



E n = -% + n, 71 = 0,3,4,5,... 



2 

where ^ n (x) is the wavefunction and E n is the energy and the polynomial P n {x) relates 
to the Hermite polynomial as follows 

1 for n = 0, 

P,,U) [ (3) 
H n {x) + AnH n ^2(%) + 4n(n — 3>)H n -±(x) for n = 3,4,5,... 



* d.agboola@maths.uq.edu.au 



1 



Also, using the supersymmetric approach [2] and for certain values of the parameters g a , 
a and u, potential ([1]) has been shown to be a supersymmetric partner of the harmonic 
oscillator potential. Moreover, in a very recent work [3], by employing the Mobius trans- 
formation, the Schrodinger equation for potential (pQ) was transformed into a confluent 
Heun equation and a simple and efficient algorithm to solve the problem numerically ir- 
respective of the the values of the parameters was presented. In addition, the 3D case 
of the potential was studied for the quasi-polynomial solutions in cases where the poten- 
tial parameters satisfy certain conditions and using the asymptotic iterative method, the 
authors obtain numerical solutions to the problem for a more general case [3]. 

On the other hand, due to the restricted number of exactly solvable systems, recent 
attentions have been on the systems with partially solvable spectral. Such systems are said 
to be Quasi- Exactly Solvable (QES). Thus a quantum mechanical system is called quasi- 
exactly solvable, if only a finite number of eigenvalues and corresponding eigenvectors can 
be obtained exactly [5]. An essential feature of a QES system is that having seperated 
the asymptotic behaviours of the system, one gets an equation for the part which can be 
expanded as a power series of the basic variable. This equation unlike an exactly solvable 
equation with two-step recursive relations, possesses a three-step recursive relations for 
the coefficient of the power series. The complexity of the recursive relations does not 
allow one to guarantee the square integrability property of the wavefunction. However, 
by chosing a polynomial wavefunction, one can terminate the series at a certain order and 
then impose a sufficient condition for the normalization. By so doing, exact solutions to 
the original problem can be obtained, but only for certain energies and for special values 
of the parameters of the problem. 

Solutions to QES systems have mostly been discussed in terms of the the recursion 
relations of the power series coefficients, which is mostly expressed in terms of the (gen- 
eralized) Heun differential equations. Although the solutions obtained in connection with 
the Heun equations are exact but the procedures involved are quite ambiguous, thus ex- 
punging the closed form of the solutions. In a series of recent studies [6]-|10j. the Bethe 
Ansatz Method (BAM) has been used in obtaining the solutions to QES systems. This 
method did not only yield exact solutions, it also preserve the closed form representation 
of the solutions. For instance, the BAM has been used to obtain the solutions of QES 
difference equation [6] and the exact polynomial solutions of general quantum non-linear 
optical models [8] and recently, the method has also been used to obtain the exact 
solutions for a family of spin-boson systems [TO]. 

The purpose of this work is to extend the study of the Generalized Isotonic Oscillator 
(GIO) to a relativistic case, with in the framework of the pseudospin symmetry Dirac 
equation, using the Bethe Ansatz Method (BAM). In section 2, we reduce the Dirac 
equation with the GIO to a QES equation. A brief discussion of the Bethe ansatz method 
is given in section 3 followed by the solution to the reduced equation. Section 4 presents 
the Lie algebraic structure of the system and then some concluding remarks are given in 
section 5. 

2 Dirac equation with the GIO 

The Dirac equation for a single-nucleon with mass \x moving in a spherically symmetric 
attractive scalar S(r) and replusive vector V(r) GIO, with h = c = 1 is written as |llj-[13j 

3 

H^(r) = £ n *(r) where H = ^ ajpj + + S(r)] + V(r) (4) 
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and E n is the relativistic energy, {ctj} and (5 are Dirac matrices defined as 

where frj is the Pauli's 2x2 matrices and is a 2 x 2 unit matrix, which satisfy anti- 
commutation relations 

aja fc + a k aj = 25jf,l 

ajf3 + f3aj = (6) 
a 2 = (3 2 ' = 1 

and pj is the three momentum which can be written as 

d 

Pj = -ify = 1 < i < 3. 

The orbital angular momentum operators Ljk, the spinor opertaors Sjk and the total 
angular momentum operators Jjk can be defined as follows: 

Ljk = ~Ljk = ^ x i~Q^, ^ Xk ~dx~' 1 = ~^ k i = ^j^k/2, Jjk = Ljk + Sjk- 

3 3 3 

L 2 = Y,L%, S 2 = ^S] k , J 2 = Y. J % Hi<K3. (7) 

j<k j<k j<k 

For a spherically symmetric potential, total angular momentum operator and the spin- 
orbit operator K = — /3(J 2 — L 2 — S 2 + 1/2) commutate with the Dirac Hamiltonian. For a 
given total angular momentum j, the eigenvalues of K are k = ±(j + 1/2); k = —(J + 1/2) 
for aligned spin j = £ + ^ and k = (j + 1/2) for unaligned spin j = £ — i. Moreover, the 
spin-orbital quantum number k is related to the orbital angular quantum number £ and 
the pseudo-orbital angular quantum number £ = £+1 by the expressions k(k+1) = £(£+1) 
and k(k — 1) = £(£ + 1) respectively for k = ±1, ±2, .... The spinor wave functions can be 
classified according to the radial quantum number n and the spin-orbital quantum number 
k and can be written using the Dirac-Pauli representation 



(8) 




where F(r) and G(r) are the radial wave function of the upper- and the lower-spinor 
components respectively, Yj m (9,(j)) and Yj m (9,(j)) are the spherical harmonic functions 
coupled with the total angular momentum j. The orbital and the pseudo-orbital angular 
momentum quantum numbers for spin symmetry £ and and pseudospin symmetry t refer 
to the upper- and lower-component respectively. 

Substituting Eq. (J8]) into Eq. ([8]), and seperating the variables we obtain the following 
coupled radial Dirac equation for the spinor components: 

<l + l) F (r) = [n + E n -A(r)]G(r) 



dr r 

±.-£)G\r)-- /,-/+,+-(r)]F(r) (9) 
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where E(r) = V(r) + S(r) and A(r) = V(r) — S(r). Using Eq. (9a) as the upper component 
and substituting into Eq. (9b), we obtain the following second order differential equations 



d 2 


k(k + 1) 


dr 2 




' d 2 


k{k — 1) 


dr 2 





[n + E n -A{r)]\}i-E n + V{r)] + 



[fx + E n -A{r)][ f x-E n + E(r)] 



dr \ dr t ) 

[fi + E n - A(r); 

( _d_ 

dr 



F(r) = 0, 



t_d k\ 

V dr r J 



G{r) = 0. 



(10) 



\H -E n + E(r)] 

To solve these equations, we employ the concept of pseudospin symmetry [121 US] in 
which V(r) + S(r) = C ps , C ps being the pseudospin constant. This implies dS J^ = ^ an< ^ 
hence Eq. (10b) takes a simple form 



f k(k - 1) 

| <ir 2 r 2 



[/x + - A(r)] [/i - £ n + ^ G(r 



0. 



If we take 



A(r) = + 2g 



r 2 - a 2 



(r 2 + a 2 ) 2 ' 

and introduce the dimensionless quantity z = f3 n ujr 2 , Eq. (11) becomes 



zG"{z) + l -G'{z) 



z k(k - 1) g(z- uja 2 (3 n ) 
4 ~Tz 2 (z + ua 2 (3 n ) 2 



G(z) = 0. 



where 



~~ fin — M ~~ -^n + Cp 



and a 



(/3 2 + 2 M + C ps ) ft, 



, ps „ 

From the asymptotic behaviour of Eq. (|13p , one may seek a solution of the form 

G(z) = (z + fi n oja 2 ) b+1 W 2 e~ z ' 2 ^(z), 

to obtain 



(11) 
(12) 

(13) 

(14) 
(15) 



zG"{z) + 



+ 2b + k 



2(6 + l)/3 n ua 2 
z + fi n uia 2 



P 2 n g/2 -(b + l)(K + b + Pnuja 2 + 1/2) 

z + fi n U)a 2 

where we have assumed 

6(6 + l)-/3 2 5 = 



G(z) 



G'(z) 



k 5 
a+6+ 2 + 4 



G{z). 



(16) 



-1- 



(17) 
(18) 



It can be checked that a further transformation 

G{z) = z v (p(t), t = z + (3 n u;a 2 

does not change the structure of the differential equation (|16p provided v = and f = 
— k + 1/2. This indicates that the solutions are of double algebraic sectors, with the 
even solution correponding to v = (k > 0) and the odd solution corresponding to 
v = — k + 1/2 (k < 0). Thus if we use the change in variable t = z + (3 n Lua 2 , Eq. (fTTj) takes 
the form 

" '5 



t(t-^ n o;a 2 )^"(t) + 



/ k 5 

t[a + b + u + - + - 
V 2 4 



+ 26 + k + 2i> + fi n ua 2 \t-t- 2f3 n uja 2 {b + 1) 



ip(t) = [fi 2 n g/2 -(b+l)(K + b + 2u + /3 n ua 2 + 1/2)] <p(t). 



(19) 
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3 The Bethe ansatz solutions to relativistic GIO 



In this section, we give a brief description of the BAM of solving QES equation and then 
use the results to solve Eq. (|19p . For interested reader, detailed account of the method 
can be found in Ref.[lQ], We consider the differential equation of the form 



S(t) = 0, 



(20) 



where P(t),Q(t) are polynomials of degree 2 and R(t) is a polynomial of degree 1, which 
we write as 

2 2 1 

P(t) = Y,Pkt\ W = R{t) = Y,r k t k . (21) 

fc=o fc=o fc=o 

where pk, qk an d r/% are constants. This equation is quasi-exactly solvable for certain 
values of its parameters, and exact solutions are given by degree n polynomials in t with 
n being non-negative integers. In fact, this equation is a special case of the general 2nd 
order differential equations solved in [lOj by means of the BAM. Applying the results in 
[TU] . we have 

Proposition 1 Given a pair of polynomials P(t) and Q(t), then the values of the co- 
efficients tq and r\ of polynomial R(t) such that the differential equation Eq. (|20p has a 
degree n polynomial solution 



S(t) = J[(t-ti), S(t) = 1 for n = 



i=l 



with distinct roots t\, . . . , t n are given by 

-n = nq 2 , 



(22) 



(23) 



-r 



i=l 



where the roots t\,t2,--- ,t n satisfy the Bethe ansatz equations 



to t A 



Qitf + qitj + g 

P2tf + PlU + Pq 



1,2, 



, n 



(24) 



(25) 



The above equations ()23f l f|25 |) give all polynomial R(t) such that the ODE (|20p has a 
polynomial solution ([2^]) . 

It is interesting to note that in line with the recent work [16j, Eqs. (|23p . (|24p and (|25p 
satisfy the necessary and sufficient conditions for the differential equation Eq. (|20p to have 
polynomial solutions. It is easy to show that the necessary condition (2.10) of [16] reduces 
to Eq. (j23|) for 03^ = and for instance, if we consider the case n = 1, then tridiagonal 
determinant (Eq. (2.11) of |16j ) for Eq. (|20p takes the form 



-ro 



-ro - qi 







-qi ± Vij- 4 9o?2 



(26) 



where we have used the necessary condition (|23p . This result can be easily obtained by 
solving for parameter t in Bethe ansatz equation (|25p and substituting the value into 
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Eq. (|24p . Thus, Eqs. (|24p and (|25p are the sufficient conditions for differential equation 
(|20p to have a exact polynomial solution (|22p . However, it is important to note that one 
of the main tasks in the application of BAM is obtaining the roots of the n algebraic Bethe 
ansatz equations (|25p . For an arbitrary n, the equation is very difficult, if not impossible, 
to solve algebraically. However, numerical solutions to the Bethe ansatz equations have 
also been discussed in many applications |17j-|21j. 

By comparing Eqs. (fTUj) and (f2"0j) . we have P2 = 1, Pi = —fi n uja 2 , qi = — 1, qi = 
(| + 26 + re + 2v + PnUJa 2 ) , q = -2f3 n oja 2 (b + 1), n = - (a + b + v + § + f ) and r = 
- [$fo/2 - (6 + 1)(k + 6 + 2^ + /3 n wa 2 + 1/2)] . Thus by Eqs. CED, CCD and (J23D, we im- 
mediately have the energy equation 

+ + £ + + 2p + C ps ) ft, + 2<Vl + 4^ (27) 

and Eq. ([Ml) yields 

n(n + 26 + 2z^ + K + /3 n o;a 2 + 3/2)-^t a = ^ -{b+l)(K+b+2v+ f3 n uja 2 + 1/2), (28) 

i=l 

with the roots t\,ti, . . .t n satisfying the Bethe ansatz equation 

2 [2b + K + 2v + f3 n uja 2 + 5/2) - g - 2/3 n o;a 2 (fr + 1) 

for i = 1, 2, . . . , n. 

It is obvious from the energy equation (|27p that one deals with solutions with positive 
energy states. Moreover, the rhs of the energy equation (|27p remain unchange for quantum 
states (n, re) and (n — 1, re + 2) thereby signifying degeneracy of the energy levels between 
these states. This energy degeneracy does not depend on the potential parameters as it 
can be seen from the numerical energy values of the ground state and some excited states 
for the exact pseudospin case (Table 1). And we also note that for a given state, the 
energy values are inveresly proportional to parameter g. Moreover, for the ground state, 
n = 0, we have from Eq. (|28p 



B 2 a 

^--{b + l){K + b + 2u + p uja 2 + 1/2) = 0, (30) 
which yields the following condition: 

(4a; V -g)fil + 2wa 2 (l + 4re)/3 + 2re(2re + 1) = for v = (31a) 

(4wV -g)Po + 2wa 2 (5 - 4re)/3 + (4k - 6) (re - 1) = for v = -re + 1/2 (31b) 
Hence the ground state wavefunctions for even sector [y = 0) (re > 0) can be written as 



Fn(r) 



Go(r) 



k / 9 ?\b _io^_ 

~r K (r z + a z ) e 2 



/ /3 ajr 3 -(26-/3 t<Ja 2 +2)r 

m 



8f 



\ r 2 + a 2 



(32) 



with the parameters satisfying Eq. (|31ap . Similarly, the odd sector {y = — re + 1/2) (re < 0) 
solutions are 

odd / /3 wr 4 -(26-2ft-/3 wa 2 +3)r 2 + (2K-l)a 2 



F (r) 
G (r) 



~ r K (r z + a z ) e 2 



(33) 



\ r [r 2 + a 2 ^ 
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with the parameters satisfying Eq. (|31bp and (3$ is related to ground state energy and 
obtained from Eq. ([2?]) as 



^ ("+£ + !) = {Pl + 2 V + C ps) A) + 2ui^l+Agpl 



(34) 



The wave functions i ? o( r ) and Gq(t) do not have nodes and so the states described by 
them are ground states of the system. 

Similarly for n = 1, the Bethe ansatz equation (|29p becomes 



(26 + k + 2v + /3i^a 2 + 5/2)*! - t 2 - 2^uja 2 (b + 1) = 0, (35) 

which yields 

^a/2(6 + i/)(26 + 2z^ + 2k + 5) - /3iwa 2 (46 + 8k + 16^ - /3ia;a 2 + 3) + (k + 5/2) 2 . 



We substitute the roots into Eq. (|28|) and solve the resulting algebraic equation to obtain 
the following condition on the parameter 

g(3f (l6w 4 a 8 + g 2 - 8guj 2 a 4 ) + I6guja 2 (3f (4wV - g) (k + 1) 

+Pf [4uj 2 a 4 g (24k 2 + 48k + 11 - 8w 2 a 4 ) - 2 5 2 (4k 2 + 8k - 9)] 
+/3f [4 5 a;a 2 (16K 3 + 48k 2 - 38k - 9) - 16w 3 a 6 (8K + 11)] 
+/3 2 [2g (8k 4 + 32k 3 + 54k 2 + 26k + 3) - 24w 2 a 4 (8k 2 + 18k + 11)] 
-4ua 2 f3 1 (32k 3 + 84k 2 + 64k + 15) - 4k (8k 3 + 20k 2 + 16k + 3) = (36a) 

g/3f (I6a; 4 a 8 + g 2 - 8gu 2 a A ) + 16gua 2 f3f (-4kwV + 8w 2 a 4 + Kg - 2g) 

+2/3 4 [g 2 (4k 2 + 16k - 2l) + 2goj 2 a A (24k 2 - 96k + 83) - 16w 4 a 8 ] 

+/3 3 [4 5 wa 2 (-16K 3 + 96k 2 - 182k + 93) - 16a; 3 a 6 (8K + 19)] 

+/3 2 [4g (4k 4 - 32k 3 + 99k 2 - 131k + 66) + 8w 2 a 4 (-24k 2 + 102k - 111)] 

+4wa 2 /3i (32k 3 - 180k 2 + 328k - 195) + 4 (-8k 4 + 52k 3 - 122k 2 + 123k - 45) = 

(36b) 

where /3iis related to the first excited energy and its obtained from Eq. (|27p as 



4oj (v + | + ~J = (ft + 2fi + C ps ) /3i + 2oj^l + 4gf3 2 . 



(37) 



Thus the wavefunctions for the first excited state for the even sector {y = 0) (k > 0) can 
be written as 



F 1 (r) 
Gi(r) 



~ r K (r z + a z ) e 2 x 

/3iwr 5 ~(3 1 u(2b-2l3 1 Lua 2 +t\+4)r 3 -[fcuia 2 (2b-(3 1 uia 2 +t\+A)-2t\ (6+l)]r 

(r 2 + a 2 ) (r 2 + a 2 - if) 



(38) 
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with the parameters satisfying Eq. (|36ap and the root given as 
k f3\ua 2 5 1 



*i = h + 2 + r ~ ± ^2~ + 4 ± 2^ 2b ( 2b + 2k + 5 ) ~ Piua 2 (4b + 8k - fruja 2 + 3) + (k + 5/2) 2 . 
Similarly, the odd sector (z/ = — n + 1/2) (k < 0) solutions are 

odd 



F 1 (r) 
Gi(r) 



~ r [r 



+ a 2 ) 



e 2 x 



(r 2 +a 2 )[/3faj 2 r 4 +/3iw(/3iwa 2 +K-t5 , -2)r 2 +K(/3iaja 2 -ff)]-2r(b+l) \ 



(39) 



r(r 2 + a 2 ) (r 2 + a 2 - tf) / 
with the parameters satisfying Eq. (]36bp and the root given as 
q = b~ +^-+ l±\V3(b + 2)(26 - 2 K + 1) - /3!o;a 2 (46 - 8k - focua 2 + 11) + (k + 5/2) 2 . 



Table 1: Relativistic energy values of the GIO for various n and k and a special case of 
C ps = 0, g = {|, 2, 4} and \x = uj = \. 



n 


K > 


9=k 


5 = 2 


g = A 




1 


0.7779142 


0.5545547 


0.4385638 




2 


1.0709485 


0.7934184 


0.6363539 





3 


1.3018103 


1.0000000 


0.8157561 




4 


1.4942015 


1.1824728 


0.9806115 




5 


1.6601782 


1.3457688 


1.1328206 




1 


1.3018103 


1.0000000 


0.8157561 




2 


1.4942015 


1.1824728 


0.9806115 


1 


3 


1.6601782 


1.3457688 


1.1328206 




4 


1.8068160 


1.4934769 


1.2738521 




5 


1.9386413 


1.6283325 


1.4049965 




1 


1.6601782 


1.3457688 


1.1328206 




2 


1.8068160 


1.4934769 


1.2738521 


2 


3 


1.9386413 


1.6283325 


1.4049965 




4 


2.0587311 


1.7524556 


1.5274033 




5 


2.1692751 


1.8675089 


1.6420843 




1 


1.9386413 


1.6283325 


1.4049965 




2 


2.0587311 


1.7524556 


1.5274033 


3 


3 


2.1692751 


1.8675089 


1.6420843 




4 


2.2718886 


1.9748119 


1.7499209 




5 


2.3677989 


2.0754239 


1.8516754 



4 Hidden Lie Algebraic Structure 

One way to understand the QES theory is to demonstrate that the Hamiltonian can be 
expressed in terms of generator of a Lie algebra 

d . 9 d r, d n ... 

J" = — , J + =t 2 —-nt, J° = t- , 40 

df dt dt 2' v 1 
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which are differential operator realization of the n + 1 dimensional representation of the 
si (2) algebra. Moreover, if we write the basic equation (|20p in the Schrodinger form 



HS(t) = -r S(t) 



(41) 



where — ro is the eigenvalue of the Hamiltonian H, then it can easily be shown that 
if r\ = —nq2, with n being any non-negative integer, the differential operator H is an 
element of the enveloping algebra of Lie algebra sl{2) 

H = J°J + P1 J J- + g 2 J+ + (gi + n - 1) J° + (q + ^f) J~ + \ (~ + qi - l) . (42) 



Thus for Eq. (|19p . we have 

d 2 

H = t(t- p n cja 2 )-^ + 



and 



- + 2b + n + 2v + f3 n ua 2 ) t - t 2 - 2(3 n ua 2 (b + 1) 



t ( a + b + v + J + 7 
2 4 



r = (6 + 1)(k + 6 + 2z^ + /3 n wa 2 + 1/2) - f n gj2 
with s/(2) algebralization 



H = J°J° - /3 n ua z J u J- - J + + [ n + 2b + k + 2v + /3 n u;a 



,2 T 



-/3 n wa 2 ( 26 + 2 - - ) J" + - ( 2b + k + 2is + /3„a;a 2 + - + - 



(43) 
(44) 



(45) 



5 Concluding Remarks 

In this paper, we have extended the works on the GIO to a relativistic case by constructed 
the Bethe ansatz solutions to GIO, within the frame work of the relativistic Dirac equation. 
We showed that the governing equation is reducible to a QES differential equation which 
has an exact solution, provided the parameters satisfy certain constraints. Unlike previous 
non-relativistic cases the quasi-exact solvability of the equations has enabled us to use 
Proposition 1 to obtain closed form expressions for the energies and eigenfunctions. It 
is interesting to note that with the limits /3 n — > 1, k — > I + 1, the spinor component of 
the wavefunction, G{z) gives the non-relativistic wavefunction of the GIO, which is in 
agreement with previous works [H El [31 IH [HJ [22] 

Moreover, we reported the existence of degeneracies between the energy levels and 
the energy is inversely related with the potential parameter g. We also showed that the 
relativistic GIO possesses an underlying sl(2) algebraic structure , which is responsible for 
the quasi-exact solvability of this model. Let us remark, however, that the existence of a 
underlying Lie algebraic structure in a differential equation is only a sufficient condition for 
the differential equation to be quasi-exactly solvable. In fact there are more general (than 
the Lie-algebraically based) differential equations which do not possess a underlying Lie 
algebraic structure but are nevertheless quasi-exactly solvable (i.e have exact polynomial 
solutions ) [10] . Finally, it is pertinent to note that our method gives a more general closed 
form expressions for the solutions, however the determination of the roots of the Bethe 
ansatz equations for higher exicted states may be a major difficulty in the application of 
the method. 
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